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Constructing reparametrization invariant metrics 

on spaces of plane curves 

Martin Bauer, Martins Bruveris, Colin Cotter, Stephen Marsland, Peter Michor 

Abstract — Metrics on shape space are used to describe deformations that take one shape to another, and to determine a distance 
between them. We study a family of metrics on the space of curves, that includes several recently proposed metrics, for which 
the metrics are characterised by mappings into vector spaces where geodesies can be easily computed. This family consists of 
Sobolev-type Riemannian metrics of order one on the space Imm(Si,IR^) of parametrized plane curves and the quotient space 
Imm(S^,]R^)/DifT(5^) of unparametrized curves. For the space of open parametrized curves we find an explicit formula for the 
geodesic distance and show that the sectional curvatures vanish on the space of parametrized and are non-negative on the space of 
unparametrized open curves. For the metric, which is induced by the "R-transform", we provide a numerical algorithm that computes 
geodesies between unparameterised, closed curves, making use of a constrained formulation that is implemented numerically using 
the RATTLE algorithm. We illustrate the algorithm with some numerical tests that demonstrate it's efficiency and robustness. 

Index Terms — curve matching, elastic metric, geodesic shooting, reparametrization group, Riemannian shape analysis, shape space 



1 Introduction 

The mathematical analysis of shape has become the focus 
of intense research interest in recent years, not least 
because of applications in image analysis and computer 
vision, where methods based on geodesic active contours 
or 'snakes' are used for segmentation, tracking and object 
recognition Ii24l l25l . Another source of applications is 
biomedical image analysis, where the study and com- 
parison of shapes form a large part of the field of 
computational anatomy JSJUl. 

For the purposes of this paper we define shape to be 
the image of a smooth closed plane curve. A slightly 
narrower definition would be to define a shape as the 
outline of a smooth, simply connected domain in the 
plane. The latter definition would exclude objects like 
the figure eigth from the class of shapes, which we allow 
in our definition. 

A key problem in shape analysis is to define a distance 
function between shapes, whose numerical computation 
is feasible and which can be used to measure similar- 
ity between shapes or as the basis for classification of 
shapes. One way to arrive at a distance function is by 
equipping the space of shapes with a Riemannian metric. 
The Riemannian metric allows us to measure the length 
of paths and we can define the distance between two 
shapes to be the length of the shortest path connecting 
them. These shortest paths are called geodesies and the 
resulting distance function the geodesic distance induced 
by the Riemannian metric. 

A shape can be represented via a smooth closed curve 
c : S*^ — > and two curves c, d represent the same 
shape, if they have the same image, c{S^) = d{S^), or 
equivalently, if there exists a reparametrization map (p g 
DifF(5^), such that one curve is a reparametrization of 
the other one, i.e. c = d o cp. Here DifF(S'^) denotes the 



group of smooth invertible maps (p : ^ from the 
circle onto itself. 

We will work in this paper with the class of regular 
or immersed curves c, that is those, whose tangent 
vector doesn't vanish, i.e., c'{9) ^ 0. Upon the space of 
immersed curves, 

Imm(5\M^) := {c G C°°(S'\m2)| c'{9) 0} , 

acts the diffeomorphism group Diff(S'^) from the right 
via {(fi, c)i—>-coif. Using this setting we can identify a 
shape with an equivalence class [c], that is an element of 
the quotient space Imm(S'^, M'^)/ Diff (S*^). This quotient 

S = lmm{S\R^)/DiS{S^) 

is the shape space of immersed curves modulo 
reparametrizations . 

To arrive at a distance function on shape space two 
steps are necessary. First we have to define a Riemannian 
metric on the space of immersed curves. A Riemannian 
metric is an inner product measuring the length of 
infinitesimal deformations of a curve. Such deformations 
h, k are represented by vector fields along the curve c 
and the inner product, which depends on the curve, is 
denoted by Gc{h,k). The first step consists of defining 
the metric and computing its geodesies on the space of 
immersed curves. 

If the metric is invariant under the action of the 
reparametrization group DifF(S'^), then it induces a Rie- 
mannian metric on shape space S, which in turn gives 
rise to the geodesic distance function. The second step 
is to find the right representatives c, d of the equiva- 
lence classes [c] and [d], such that the geodesic distance 
dist^"™(c,d) coincides with dist'^([c], [d]). 
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1 .1 Shape metrics and related work 

The simplest reparametrization invariant metric on 
Imm(S'^, M-^) is the L-^ -metric 

Gcih,k) = I {h,k)ds, 

where we integrate over arc-length ds — \c'{6)\d9. How- 
ever it was found in ||T6|| that the geodesic distance 
induced by this metric vanishes, i.e. the distance between 
any two shapes is 0, which renders this metric unsuitable 
for shape analysis. 

A possible way to overcome this is to add terms 
involving higher derivatives of h and k to the metric, 
like 

Gc{h,k)^[ {h,k)+A{Dsh,Dsk)ds, 

where Dsh = -^^h' denotes the arc-length derivative 
of h and A > is a constant. This leads to the class 
of Sobolev-type metrics introduced in 1191 and studied 
further in |3, 21 J. 

Another family of metrics, the almost local metrics 
ISI ITtII , prevent the geodesic distance from vanishing by 
introducing a weight function in the integral. Examples 
of weight functions involving the curvature or length are 
w{e) = 1 + Ak{9)^ and w{e) = £{c)i. 

Sobolev-t5rpe metrics of arbitrary order were studied 
in ||5l [15j |T9| . Although they are a natural generaliza- 
tion of the L^-metric from a theoretical point of view, 
their numerical treatment is rather involved. This stems 
mainly from the fact that the geodesic equation of a 
Sobolev-type metric of order fc is a highly nonlinear 
PDE of order 2k. While there have been some attempts 
to solve the geodesic equation directly for order 1 for 
curves |20J and surfaces |2J, the metrics of higher order 
are still practically untouched. 

To bypass these difficulties, one can restrict the atten- 
tion to Sobolev-type metrics of order one, in particular 
to the family of metrics studied in 1,20 J , which are of the 
form 

GJ?'^(/i,/i) = / a^{Dsh,nf + b^{Dsh,vfds, 

with constants a,b E M+ and with v and n denoting the 
unit tangent and normal vectors to c. Metrics of these 
form are called elastic metrics. The term involving the 
normal vector can be seen as measuring the bending of 
the curve c under the deformation h, while the other 
term measures the stretching of c. 

Following ideas of |26|, it was shown in [27J that it 
is possible to find explicit formulas for geodesies of the 
elastic metric, when = a^. To achieve this a curve c 
was represented by the square-root of its velocity vector, 
y/d, with c' being interpreted as a complex number 
In this representation the G^' ^-metric has a particularly 
simple form which allowed to find explicit formulas for 
geodesies between two curves as well as the geodesic 
distance. 



A similar motivation underlies the introduction of the 
i?-map in [23 J. The i?-map is a transformation, which 
maps a curve c to R{c) = c! j \/\d\ with the effect that the 
G" '' -metric simplifies for 46^ ~ a} . While not allowing 
to write down explicit formulas for geodesies, this map 
greatly simplifies their numerical computation. 

Related to the i?-map is the Q-map, introduced in IITTl , 
which maps a curve c to (5(c) ~ yj\d\c. The advantage of 
this map is that it generalizes easily to surfaces. It does 
however present theoretical difficulties as explained in 
Section |5] 

There are other approaches to define a Riemannian 
metric on shape space: the Large Deformation Diffeo- 
morphic Metric Mapping approach fS", '6], where a Rie- 
mannian metric is induced from the diffeomorphism 
group DifF(]R^) of the ambient space or use of conformal 
welding f22] to represent shapes as diffeomorphisms of 
the circle. 

1.2 The reparametrization group Diff(S'^) 

After equipping Imni(5^,M^) with a Riemannian metric, 
there are various possibilities how to perform the mini- 
mization over the reparametrization group. The problem 
is challenging because of the nonlinear nature of both 
DifF(S'^) and since all the spaces involved are infinite 
dimensional. 

One approach is to replace the infinite dimensional 
space DifF(S'^) by a finite-dimensional one and to per- 
form the minimization over this smaller space. A pos- 
sible choice for the smaller space, used in fSS^, are 
diffeomorphisms 93, for whom the Fourier series of yV' 
is truncated at a fixed length. 

In |6| elements of DifF(S'^) were generated as flows 
of vector fields. The advantage of using vector fields 
is that they form a linear space, which enables us to 
use gradient-based optimization algorithms. However, 
due to the regularization terms needed to ensure con- 
vergence of the algorithm, the computed distance failed 
to be symmetric. 

In the paper |24] the authors describe an iterative 
procedure directly on Diff(S'^) based on geometric con- 
siderations. We will expand on this minimization scheme 
in Section |6] 

1 .3 Overview of the paper 

In this paper we introduce a i?-transform, which enables 
us to represent all elastic metrics satisfying 46^ > as 
pullbacks of the flat L^-metric, generalizing the result of 
|23], which only applied to 4fe^ — a? . 

Then we proceed to analyse the mathematical prop- 
erties of this metric, first on open curves in Section |2j 
followed by closed curves in Section |3] In particular we 
derive a formula for the geodesic distance between open 
curves and show that the space of open curves is flat in 
the sense of Riemannian geometry. We show in Section 
|4]how the latter result implies that the sectional curva- 
ture on shape space of open curves is non-negative. In 
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addition we derive a formula for the sectional curvature 
on the space of closed parametrized curves. 

Section |5] is devoted to general shape transformations, 
that can be used to construct reparametrization invariant 
metrics on Imm(S'^,M^). Both the Q-map from flT] and 
the representation of curves used in l,27J fit into this 
framework. 

We move on in Section |6] from theoretical considera- 
tions to numerical computations. First we propose, using 
the i?-transform developed in Section |2| an efficient 
shooting algorithm for parametrized curves. Then we 
show how the method of optimizing over reparametriza- 
tions, that was also considered in 1241 , can be seen as a 
gradient descent with respect to a certain Riemannian 
metric on Diff(S'i). 

Finally in Section [7| we show some experiments and 
identify the limits of the metric for practical use. 

2 The i?-TRANSFORM FOR OPEN CURVES 

We want to study a family of Riemannian metrics on the 
space of open and closed curves, which are constructed 
as pullback metrics of a simpler metric under a transfor- 
mation of the curve. Let us define the i?-transform of a 
plane curve c e lmm([0, 27r],M^) by 

R"-^ : Imni([0,27r],M2) ^ C°°([0, 27r], K^) 



Here a, 5 e with 46^ > a} are positive numbers, 
whose meaning will become clear, when we compute the 
pullback metric. We will omit the parameters a, h when 
the meaning is unambiguous. We denote by w = the 
unit length tangent vector of c. The i?-transform maps 
an open plane curve to a space curve. We will see that 
equipping the space C°° ( [0, 27r] , W") with a flat L^-metric, 
i.e. considering the vector space C°°([0, 27r],M'^) with 
the i^-rnner product as a Riemannian manifold, will 
generate Sobolev metrics of order one via the pullback by 
the i?-transform. For the choice of parameters 4&^ = d} 
the i?-transform reduces to the i?-map as studied in fSS]. 
The metric induced by the i?"' ''-transform also coincides 
with the metric studied in EOll corresponding to the 
parameters ^/a, \/h. 

First let us note the following properties of the R- 
transform. 



Lemma 2.1. For a curve c e lmm([0, 27r], M^) 
dijfeomorphism ip e Diff^([0, 27r]) we have: 

• Equivariance under reparametrizations 

R{coip) = \ip'\^^^.{R{c) o ip) 

• Translation invariance 

R{c + p)^ R{c) for peR^ 

• Scaling property 

R{p.c) = p^/^.R{c) for p £ M>o 



and a 



Preservation of length under reparametrizations 



\R{cop)Yde 



SI 



\R{c)Yd6 . 



Proof: We will show the first property, since it will 
be important for the reparametrization invariance of the 
pullback metric. From (c o p)' = [d o p)p' it follows 
that the unit length tangent vector v behaves under 
reparametrizations as 



v{c o ip) 



{cop)' 

\{cop)' 



p — O = V{C) O p> 



The term |c'|"'"/^ reparametrizes as |(co(^)'|i/^ — \p'\^/'^\c'o 
p\^^'^, and we have established the first property. The 
other properties can be verified by a simple calculation. 

□ 

In the next theorem we calculate the pullback of the 
flat i^-metric under the i?-transform. 

Theorem 2.2 (The pullback metric on hxim). The pullback 
of the L'^-inner product on C°°([0, 27r],IR'^) to the manifold 
of immersions by the R''''^ -transform is given by: 



Jq 



DM, v)"^ ds. 



where n is the unit length normal vector to the curve c, the 
arc-length derivative of h along the curve c is denoted by 
Dgh = and ds = \c'\d9 is the integration with respect 
to arc-length. The metric can be written in the form 



G:'\h,h) 



{P^'''{h),h) ds 



where the associated differential operator P^''' is given by: 

P^'\h) =-a'^{Dlh,n)n-h^{Dlh,v)v 

+ (a^ - h^)K{{D,h, v)n + {D,h, n)v) 

+ (<527r - 5Q){a^{n,Dsh)n + b^{v,Dsh)v), 

with K being the curvature of the curve. 

The pullback defines a Riemannian metric on the space 
lnim([0, 27r],M^)/ trans of plane curves modulo translations. 

The metric G°-'^ is invariant with respect to the 
reparametrization group Diff^([0, 27r]). 

Proof: The pullback metric is defined via 

Gt-\Kh) = {T,R{h).T,R{h))^,, 

and hence we need to compute the derivative of the 
i?-transform. First we compute the derivatives of the 
functions c ^ |c'|^^^ c H> which are 



Tcv{h) = ^ADsh,v)c! + ^h' 
\c'\ \c'\ 

— Dgh — {Dgh, v)v — {Dgh, n)n . 
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Now computing the derivative of the i?-transform is 
simple. 



r.'|l/2 



a{Dsh,n)n + ^{Dsh,v)v 



lVW^{D,h,v) 
Therefore the puUback metric is 



G^/{h,h)= r 

Jo 



a{Dgh,n)n + -{Dsh,v)v 



+ 



+ {b'-^]{D,h,vy]\c'\de 



/•2-ir 

= 

Jo 



{Dgh,nf + b'^{Dsh,vf ds . 



To obtain the formula for the operator P^-^ we apply 
integration by parts to the above formiila: 

/ o^iDsh, {Dsh, n)n) + b^{Dgh, {D^h, v)v) ds 
Jo 

= {a^{h, {D,h, n)n) + h'ih., Ds{Dsh, l^" 



Jo 



{h, D, {(D,h, n)n)) + b''(h, D, {(D,h, v}v)} ds 



{o?{h, {Dsh, n)n) + b'^{h, D^iDsh, v)v)) 



|27r 



/■27r 

/ a^(/i, (-Dg/i, n)n — {Dsh, Kv)n — {D^h, n)Kv) ds 
Jo 



-I 

Jo 



b^{h, {Dgh, v)v + {Dsh, Kn)v + {Dsh, v)Kn) ds . 



Here we used that k = {DsV,n) = —{Dsn,v). 

The metric G"^'^ is invariant under reparametrizations, 

G'i;^^{hoip,hoip) = G-/{h,h), 

since it it written in terms of operations Ds,v,n which 
are equivariant with respect to reparametrizations. □ 

2.3 Image of the P-transform 

To characterize the image of the i?-transform we note 
that we can reconstruct c' from the first two components 

of R{c) via 

and hence a|c'|^/^ = \/R\{c) + i?2(c). This implies the 
relation 

which can be written in the form 

(462 - a^) [R^{cf + R2{cf) = a^Rsicf . (1) 
Let us define the following cone 
C^'" = {q€M.^ : (462 - a^){ql + g|) = a^ql q, > o} 



in R^. Then the image of the i?-transform consists of 
curves, which lie in C"'^. That is for a curve q e 
C°°([0,27r],M3) we have 

qGimR-^ q(e) e C"■^ G [0, 2n] , 
imi? = C°°([0,27r],C"'''). 

In the special case 46^ — the i?^"''' -transform has no 
third component and the image of the i?2a,a-t]-ansform is 
the open set consisting of all curves, which avoid (0, 0) e 
IR2. The latter condition arises from the requirement that 
the curves be non-degenerate. 

The inverse of the i?-transform can be computed using 
the identity 

" 2a6' ' {R2{c) 

Therefore 

Ri:imR^ lmm([0, 27r], K^)/ trans 



Ri{q){0)=po 



1 r 

2abJo 



q2{e) 



d0. 



Of course the inverse of the i?-transform is defined only 
up to translation, which manifests itself as the freedom 
to choose the starting point po £ of the integration. 

The cone C"'^ is a flat hypersurface in R"^, since there 
is a isometric covering map from the polar coordinate 
domain {{r,ip) : r > 0} to C"'^ given by 

'^cos(m(^)\ 
Q{r,'P) = ( ^ sin(my ) , 



2b 



where we introduced the constant m = ^. A straight- 
forward computation 



{4l,'-a')iql+q^) = 



462 - a2 



2 2 

-r =a q^ 



checks that the image of this map is the cone C"'*'. From 



dq? + dql + dql = ^dr^ + r'^d(f'^ + 



46^ - a' 
462 



dr' 



dr^ + r^dip^ 



we see that the map q{r, ip) is an isometry from the 
Euclidean metric in ¥>?, which has the expression dr"^ + 



r'^dip^ in polar coordinates, to the natural metric on the 
cone C"'^ 

The inverse map is determined only up to a multiple 
of 27r and is given by 



r{q) = 

¥'(a) 



26 



?93 



7462 

— ( arctan ( — 
26 V \qi 



2fc7r 
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with k E Z. Using the inverse map we can write the 
distance function on the cone 

dist\q,q) = {xi{q) - xM? + {^^{q) " X2{q)f (2) 
= r{qf + r{qf - 2r{q)r{q) cos {ip{q) - ip{q)) 

= mm — ^ \ Qr, + qr, ~ 

— 2(7c.(7„ cos I — I arctan — — arctan — ) H — kn 
' \2b\ qi qJ b 

The minimum appears, because the angle is only deter- 
mined up to a multiple of 2tt. 

Theorem 2.4. The metric G"-'^ on open curves is flat. 
Geodesies are the preimages under the R-transform of geo- 
desies on the flat space imi?. 

A path of curves q : M x [0, 27r] — > G"-'^ in im R is a 
geodesic, if for each 9 e [0, 27r] the curve t q{t,9) is a 
geodesic in C'^. 

The geodesic distance between c,c E lmm([0, 27r],Ili2) is 
given by the integral over the pointwise distance, 

dist{Ric){e),R{c){9))dd. 



However the minimum over k E Z is not to be taken 
pointwise, but only once for all values of 9. This corresponds 
to choosing a continuous lift of the curve R{c) via $. 

Proof: Since the cone C" '' is flat in the sense of 
Riemannian geometry, so is the space C°°([0, 27r], C"'*") 
of curves in C" '' with respect to the L^-metric given by 

Gg{h,k) = / {h,k)d9, 
Jo 

for q e C°°([0,27r],C"'''') and h,k tangent vectors at q. 
Note that this metric does not depend on the basepoint 
q. The metric is the same at all points in the space. Note 
also that the image of the i?-transform 

imi? = C°°([0,27r],C°'^), 

equals the set of curves, which lie in the cone C"'^. It 
is a property of the L^-metric that geodesies in im R are 
given by paths of curves q{t, 9), such that for each fixed 
9, the curve q{-,9) is a geodesic in the cone C" ''. The 
length of this geodesic will be given by an expression of 
the form Q with some fc e Z, not necessarily the smallest 
one. The length of the path q{t, 9) in im R is given by the 
integral over the lengths of each point-wise path q{-,9). 
Since we have a continuous family of geodesies on C"'^, 
the value for k will be the same for all 9 £ [0, 27r]. Hence 
the geodesic distance between two elements q,qE im i? 
is given by the integral 



2tt 



dist{q{9),q{9)d9 , 



where the minimum over fc e Z is taken only once for 
all values of 9. □ 
See \7j Theorem 9.1] for more details on the flat L^- 
metric. 



3 The i?-TRANSFORM FOR CLOSED CURVES 

In this section we want to consider the i?-transform 



acting on closed curves. First note that Lemma 2.1 and 
Theorem 2.2 remain valid, if we replace open curves by 
closed ones. 

Restricting our attention to closed curves C°°{S^ 
means that we impose additional constraints on the 
image of the i?-transform. From the inversion formula 



Rl{q){9) - 



1 



we see that a curve q is the image of a closed curve only 
if the condition 



Jo 



q2{0) 



d9 = 



is satisfied. Let us denote the image of the i?-transform, 
restricted to closed curves, by 

^^'b = {q(z C°°(S'\ C"''') : Fiq) = 0} . 

Note the difference between the cone C" '', which is 
a submanifold of M'^ and the space '^"■'b, which is a 
submanifold in the space of curves. 

Theorem 3.1. The image '^"■^^ of the manifold of closed curves 
under the R-transform is a codimension 2 submanifold of the 
flat space C°°(S'\ C'^^''). 

A basis of the orthogonal complement (Tg*^"'^)^ is given 
by the two vectors 



Ui{q) = 



U2{q) = 




v/462 - a2 



v/452 - a2 



0^ 

\12J 



Proof: A basis of {Tq'^if"-'')-^ can be computed by pro- 
jecting the gradients of two components of the function 
F = (Fi,F2) to the tangent space of C°°(5\ C"^-''). Let 
q e "^"'^ be a curve and h e C°°(S'^, M^) a tangent vector. 
Then 



27r 



glfel + 92^2 

\/ql + ql 

q\h\ + q^hi 



Thus the two gradients are 

grad^Vi(g) = 

grad^V2((7) = 



92 



q\ + q| h-Y d9 



12 



ql + ql h2 d9. 
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Differentiating the governing equation for the cone Q 
we obtain that the tangent space TqC°°{S^,C'^''') is given 
by all curves h, which satisfy the pouitwise condition 



(46^ - a'^){qihi + (72/12) = a^qs^s ■ 



(3) 



A projection (not neccessarily orthogonal) of the vectors 
can be found by choosing such that the 
vector grad^ Fiil) + (0,0,1)^ satisfies l|3|. A simple 
computation shows that for the component Ei the right 
value is \i — 2qi ^^b^-a ^^ jj^jg completes the proof. □ 
We can now use the Gram-Schmidt procedure to 
compute an orthonormal base of (T^ 'if '')-'-, which is 
better suited for computational purposes. The formulas 
however do not reveal more structure. 

3.2 The curvature of the space ° 

As mentioned previously the i?-transform has no third 
component in the case b = 2a and the cone C^" " reduces 
to the space \ q. Let Ui{q),U2{q) G C°°{S\R'^) be the 
orthonormal basis of (TgC^" ")-'- derived via the Gram- 
Schmidt procedure from the basis Ui{q),U2{q)- Using 
this ONB we can express the curvature of the cone ^. 

To do so we take a constant vector field q {q, h) on 
C°°(5'^,M^) and its orthonormal projection 

X,M) = h- {Ui{q),h}U[{q) - {UM,h)UM e Tq^^'-''^ . 

Then we take the flat co variant derivative in Imm(S'^ , M?) 

Vx,M)Xk{q) = Dq,x,^q)Xk{q) = d{Xk{q)){Xh{q)) . 

The orthonormal projection of this vector field onto 
{Tq'^'^'^'"-)-^ is then equal to the value of the second 
fundamental form S e (r^'r^'^'")-^, i.e.: 

S{Xh{q),Xk{q)) = {d{Xk{q)){X^,{q)),^J,{q))L^Ul(q) 

+ {d{Xk{q)){Xn{q)),U2{q))L-U2{q) ■ 

The curvature of "^^a.a ^ jg then given by the Gaufi- 
equation UHl 26.4]: 

{R{Xh (q) , Xk {q))Xk (q) , X,, (q)) = 

= -\\S{X,,{q),Xkiqmh 

+ {S{Xk{q),Xk{q)),S{Xj,{q),X^{q))}L2 

= -{d{X,M){Xk{q))Jjl{q))h 

- {d{Xk{q)){Xn{q)),U2{q))l^ _ 

+ {d{Xn{M^hm.VM)L-{d{XM){Xk{q)),UM)L- 

+ {d{Xhiq)){Xh{q)),U2{q))L-^{d{Xk{q))iXk{q)),U2{q))L^. 

4 The INDUCED METRIC ON SHAPE SPACE 5 

In the rest of the paper we will mainly consider closed 
curves. However the results can be easily reformulated 
for the case of open curves. 

The shape space S denotes the space of un- 
parametrized plane curves, which can be represented as 
the quotient S := Imni(S'^,M^)/Diff(5^) of parametrized 



curves modulo parametrizations. Associated to this quo- 
tient is the natural projection 

TT : hTini(S'\M2) ^ Imm(S'\M2)/DifF(S'i), 



which maps a curve c to its image C 

Qa.b 



By Theorem 



2.2 



7r(c). 
on Imm(5^,I 



the metric U"'" on lmm^^",JK'j is 
invariant with respect of the reparametrization group 
DifF(S'^). Therefore there exists a unique Riemannian 
metric G on shape space, such that the projection tt is a 
Riemannian submersion. Associated to the projection tt 
the decomposition of the tangent bundle rimm(S'^,M^) 
into horizontal and vertical parts. The vertical bundle 
Ver is the kernel of the projection tt, i.e. Vcr = kerTTr, 
and the horizontal bimdle 

Hor(c) = Ver(c)^ C TJmni(S'\ M^) 

is defined as the orthogonal complement of Ver with 
respect to the Riemannian metric G" ''. The action of 
DifF(S'^) on Imni(S'^,M^) induces an infinitesimal action 
of its Lie algebra X{S^), given by 

with /i e This defines a vector field on the 

space Imm(S'^,M^). The vertical bundle consists of the 
image of all infinitesimal vector fields 

Ver(c) = {C^(c) : ^i e X{S^)) . 

When we fix a curve c, there is a one-to-one correspon- 
dence between the vector fields on the circle X{S^) and 
the space Ver(c) given by the map i-> Cp(c). 

From the theory of Riemannian submersions [18] it 
follows that 

• Geodesies on shape space with respect to G cor- 
respond to horizontal geodesies on the manifold 
Imni(5'^, M^) of parametrized curves with respect to 
QaM^ Horizontal geodesies on Imni(S'^ , M^) are those, 
whose tangent vector lies in the horizontal bundle, 
i.e. dtc{t) e Hor(c(i)). 

• The geodesic distance on shape space can be com- 
puted using the formula 

distfCo, Ci) = inf dist(co, ci o 05) 
¥3eDiff(si) 

where dist on the right hand side denotes the 
geodesic distance on Inim(S'^, M^) with the G"'^- 
metric. 

• The curvature of the shape space can be calculated 
using O'Neil's curvature formula, see for example 
IflSil . Given two orthonormal vector fields X, Y on 
the space S of unparametrized curves, the sectional 
curvature K is given by 



KsiX,Y) 



,(x,r) + ^|[x,rr'-f . 



Here X, Y are horizontal lifts of the vector fields 
X, Y to Imm(S'\ M^) and [X, F]™"^* deiiotes the ver- 
tical projection of the vector field [X, F]. 
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Theorem 4.1. Consider the shape space of open curves S = 
lmm([0, 27r],IR^)/DifF([0, 2tt]) with the metric that is induces 
from the G^'^-metric on parametrized curves. The curvature 
of this space is non-negative. 



Proof: Theorem |2.4| shows that the space 
lmm([0, 27r],M^) of parametrized curves is flat, which 
implies Kii^a\{X,Y) = in O'Neil's curvature formula. 
The only remaining term 



Ks{X,Y) 
is clearly non-negative. 



vert 1 2 



□ 



5 Relating shape transformations 

In this section we describe a general method of construct- 
ing reparametrisation invariant metrics on the space 
of plane curves. This method will include the metrics 
studied by Srivastava et al. fTS\, Mio et al. |20|, the 
analogue for curves of the surface metric from Kurtek 
et al. llOl [ni and the method of Younes et al. [27] as 
special cases. 

Let us consider a general transform F : Inim(S'^, M^) ^■ 
C°°(5^,M") mapping plane curves to curves in M" for 
some ri e N. The transform F induces a reparamet- 
risation invariant metric, if it satisfies an equivariance 
property. 

Lemma 5.1. If F : Inim(S'i, M^) ^ C°°(5\M") satisfies 
the equivariance property 



F{c o ip) — ip'^F{c) o If, 



(4) 



with c e Imm(S'\M^) and (p € T>\&{S^), then the pullback 
by F of the flat L'^-metric on C°°(S'^,M") fs invariant under 
the reparametrization group Diff(S'^). 

Proof. Let us denote the pullback metric by . It is 
given via the formula 



G',{h,h) = / \T,F{h)Yd6. 

The infinitesimal version of the equivariance property is 

T,ovF{hop,) = ^'^T,F{h)o^ 
and therefore we see that 

Gi^{hocp,hoip)= f \T,o^F(ho^)\^de 

\TrF{h)o^\''\^'\de = G^AKh), 

i.e. the pullback is invariant under DifF(S'^). □ 
If additionally the function F is infinitesimally injec- 
tive, i.e. TcF is injective for all c e Imm(S'^, M^), then 
the two-form G^ is positive definite and F thus induces 
a reparametrisation invariant Riemannian metric on the 
manifold of immersed curves. 



5.2 The Q-transform 

Another member in this family of shape transformations 
is the Q-transform, which was introduced in Mani et al. 
IIT4II for curves and in Kurtek et al. IITOl [TTI for surfaces. 
It has properties similar to the i?-transform, and yet is 
much more difficult to study. 

Following IITH let us define for a curve c e 
Inim(S'^, R-^) the Q-transform by 

Q : Imm(S'\]R2) ^ C°°(5\E2) 
Q{c)^\c'\h. 

First let us note the following properties of the Q- 
transform. 

Lemma 5.3. For a curve c G Imm(S'^, M^) and a diffeomor- 
phism ip G DiflF(S'^) we have: 

m Equivariance under reparametrizations 

Q(co^) = |^f/2.(Q(c)o^) 



• Behavior under translations 

Qic + p) =Q{c) + \c'\^^^pfor peR'' 
m Scaling property 

Qip.c) = p^l^.Q{c) for p e M>o 

* Equivariance under rotations 

Q{e'"c) = e'^Qic) for a e . 

Proof: This lemma can be proven by direct calcula- 
tions similar to Lemma l2.ll □ 
Note that the behavior of the Q-transform under trans- 
lation of the curve is rather unintuitive. The curve is 
being translated by a vector, which depends on the local 
length of the curve. 

Theorem 5.4. The Q-transform induces a reparametrization 
invariant metric on Imm(S'^,K^), given by 

Gc{h,k)^ [ {h+ UD,h,v)c,k+ UDsk,v)c)ds . (5) 

Proof: We need to show that the metric is non- 
degenerate, i.e. that Gc{h,h) ~ only when h = 0. 
This condition is equivalent to the transform Q being 
infinitesimally injective, i.e. TcQ.h = only when h = 0. 
In the proof of Theorem 2.2 we computed 

r,(|cf/2) {h) = l\c'\'/'{D^h,v) 

and therefore we have 

T,Q.h^\c'\'/^h+^{Dsh,v)c) . (6) 

Since we are only considering immersions for which 
|c'| 7^ 0, the condition TcQ.h = is equivalent to 

{h',c')c + 2\c'fh = 0. (7) 

For all 9, where c{6) — this already implies h{9) — 0. On 
the open set {9 : c{9) ^ 0} we expand the tangent vector 
h into a part that points along c and a part orthogonal to 
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it, h = hic + /i2C^ with c _L c^. Then Q can be rewritten 



as 



{h',cJ)c + 2|cf /iic + 2\c/\^h2C^ = . 



This implies that h2 = everywhere and we are left with 
the equation 



c,c')h[ + 3\c'\-'hi = 0. 



(8) 



As before for 9 where {c{6),c'{6)) = 0, we infer that 
hi{6) = 0. It remains to study on the open set, where 
neither c{9) nor {c{6),c'{9)) vanish. This open set is the 
union of disjoint open intervals. Let us denote one such 
interval by {9q,9i). On this interval every solution of l|8| 
is given by 

with 9 e (6*0,6*1). Since we are looking for one smooth 
solution hi on and we already know that for 6*0 and di 
the solution has to satisfy hi{9o) = /ii(^i) = 0, because 
these points lie outside the open set, we see that only the 
solution with C — can satisfy this. Therefore hi = and 
with it /i = on all of 5*^ 

The expression ||5} for the metric Gc{h,k) is an imme- 
diate consequence of l|6|. The reparametrization invari- 
ance of the metric is either a consequence of Lemma 5.1 
or can be seen directly from ||5|. □ 

The injectivity of the map TcQ is essential in order 
to use the metric for shape comparisons. However we 
could not find a proof for this anjrwhere in the literature. 

Questions that were comparably easy to answer for the 
i?-transform are much more difficult for the Q-transform. 
To our knowledge it is not known how image of the 
Q-transform on either open or closed curves looks like; 
whether it is open or a smooth submanifold of the space 
C°°(5^,M^). Even finding a numerically efficient way to 
invert the Q-transform presents difficulties. 

The main reason for these difficulties is that the Q- 
transform scales the curve c with the object |c'|^/^, which 
geometrically is the square-root of the volume-element 
on the curve and should multiply tangent vectors along 
c, not the curve itself. 



5.5 The transform of Younes et al. fTT] 

In this section we will show that the method applied 
in Younes et al. [,27] to study a Sobolev-type metric on 
curves also fits into the setup described in this paper. 
The basic mapping considered in 112711 is given by 

$ : C°°([0,27r],M2\0) ^ lmm([0, 27r], M^) 



1 



qliu)~qj{u) 
2qi{u)q2{u) 



du . 



The basic mapping has the property that it pulls back 
the Sobolev metric of order one 

Gc{h,h) := / {Dsh,Dsh)ds 

JO 



to the flat L^-metric on the space C°°([0, 27r], \ 0). To 
fit this mapping into our framework we define the R- 
transform as the inverse of $, 

R : Imni([0,27r],M2) ^ C°°([0, 27r], \ 0) 



nf 



COS ■ 



\^sin ^ ^ 

where a{9) is the turning angle of c defined via 

c'{9) = \c' {9)\ {cos a{9), sin a{9)). 

Theorem 5.6 (From ||27|). The R-tmnsform induces the 
Sobolev metric of order one 

Gc{h,h) := I {D,h,D,h)ds 
Jo 

on the space of parametrized curves. 

6 Numerical computation of geodesics 

In this section we will describe a way to numerically 
compute the shortest path in shape space of closed 
curves between two shapes Co and Ci. We will do so 
for the special case 46^ = = 1. In this case the R- 
transform maps plane curves into plane curves, 

R : lmm{[Q, 2Tr],E.^) ^ {[0,2Tr],R^) , R{c) = \c'\^/'^v. 

The cone (7'^ "/^ regarded as a subset of simplifies to 
(ja.a/2 = ^2 0. Therefore the image of the i?-transform 
of open curves is the set 

imi? = C°°([0,27r],M2\o) 

of all curves, which avoid the origin, which is an open 
subset of the space of all curves. Restricting ourselves 
to the closed curves, we obtain from Theorem |3.1| that 
the image of the _R-transform is a codimension 2 
submanifold of C^iS^R^). 

6.1 Geodesies between parametrized curves 

Geodesies in the preshape space Imm(S'^, M'^) of 
parametrized curves correspond under the i?-transform 
to geodesics on '^/f with the Riemannian metric induced 
from the flat L^.metric on C°^{S^,R'^). 



Using Theorem 3.1 we can implement a projection 
operator which, given a curve q E and a tangent vector 
p to q, computes its orthogonal projection onto T^*^. 
function PRO^{q,p) 

Ui,U2^ ONB of (T/#')-L 
return p - {p, Ui)Ui- (p, U2) U2 
end function 

The forward computation of the geodesic from a curve 
q E ^€ with initial velocity p e Tq€ can be seen as a 
constrained optimisation problem. A geodesic on is a 
minimum of 



1 



u(t)\^ + m - u{t)) + X{t) ■ F{q{t)) dt . 
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Since we use the Euclidean metric on the ambient space, 
we will not distinguish between velocities u{t) and 
momenta p{t). The variable X{t) e is a Lagrange 
multiplier, which enforces the constraint 



F{q)= [ \q{e)\q{9)de, 



defining the cone ~ Fi{Q). These are holonomic 
constraints and there are already several methods, that 
can solve this problem; see e.g. [12] for an overview. We 
have chosen to use the RATTLE algorithm [Ij, which 
is a symplectic integrator that preserves the constraints 
through an explicit projection of the velocity vector, so 
that it lies on the constraint surface. The variational 
equations in continuous time have the form 

m = pit) 

We denote by N the number of time-steps used by the 
integrator, so At = 
function Exp(g,p, iV) 

90,^0 ^ q,p 

for i ^ 0, iV - 1 do 

qi+i ^ qi + Atp 

estimate initial value for A to enforce constraint 

iteratively adapt A so that g^+i e 

Pi+i ^p + VF(q,+i) 

/I chosen such that p^+i € Tq._^^^ 
end for 
return q,p 
end function 

The next step is to solve the boundary value problem 
for parametrized shapes. Given two curves go: 91 G we 
need to find the initial velocity p, such that the endpornt 
Exp^^ p is close to qi. We do so via a fixed-point iteration, 
using the fact, that our space is a submanifold of a flat 
space. 

function LOG{qo,qi, N) 
p^ AtPRO]{qo,qi- qo) 
q^EXF{qo,p,l) 
p^ FROKqo,q-qo) 
while \EXF{qa,p,N) - qi \ > e do 
q^EXF{qo,p, N) 
p ^ p + aPRO]{qo,qi - q) 
end while 
return p 
end function 

To construct the initial guess for p we project the straight 
line between go 9i to Tq^'^. This would be a valid 
initial guess. Based on experiments, however we found 
that computing one step of the forward-shooting algo- 
rithm to obtain q and projecting the straight line q — qa 
back to Tgo*^ leads to a better initial guess. The iteration 
consists of computing the endpornt q — Exp^^^ p of the 
geodesic with with initial velocity p and updating p in 
the direction qi — q. 



The parameter a controls the step-size of the iteration. 
It can be either fixed for the whole minimization or 
chosen adaptively. In our experiments we chose a as 
large as possible, while still ensuring that the distance 
decreased, compared to the last iteration. 

It is also possible to view this iteration as an approxi- 
mation to a gradient descent algorithm on the geodesic 
distance ^ 

E{p) ^ -dist(Expg^p,gi)2. 

The derivative of this function is 



TpE{Sp) = G'exp,^ p [- LogE^p^^^ p 91, Tp Exp^^, (dp) 

If we approximate the logarithm by the straight line, 
LogExp p 91 ~ 91 ^ ^^Pgo P' differential of the expo- 
nential map by the identity, Tp Exp^^^ ((5p) « Sp and the 
parallel transport along the geodesic 1 1-> Exp^^^ {tp) from 
Expg^^ p to qo by the projection to T^^, then we obtain 

TpE{Sp) w Gq„ (- Projg^ {qi - Exp^^ p) , dp) 

and hence the approximation of the gradient is given by 

VpE « - Proj^„ (gi - Exp,^ p) . 



6.2 Finding the optimal parametrization 

As explained in Section |4j finding geodesies between un- 
parametrized shapes corresponds to finding horizontal 
geodesies or equivalently finding the minimum of 

dist(C. D)= inf distfc, d o -0), 

^eDifF(Si) 

with C = 7i'(c), D ~ 7r{d) and where dist on the right 
denotes the geodesic distance on parametrized curves, 
computed as described in Section 6.1 We will compute 
this minimum using a gradient descent algorithm on the 
function 

E{ip) = -dist{c,doipif . (9) 



We use doLpi instead of docp, in order to have a left action 
of Diff(S'^) on the space of curves. We want to compute 
the right-trivialized gradient W^pE of E, defined as 

(V^£:, M)3e(Si) = T^E{fi. o (f) 

for /i e X{S^) and for a choice of an inner product 
on X(5i). 

From Riemannian geometry we know that the gradi- 
ent of 



is given by 



G{p) = idist(p,<7)2 



VpG = - Logp q . 



Hence 

T^E{S(p) ^ Gdo^ii" Log^„^, c,-{do ipi)' Sip o ipi) 
and by writing 5ip ~ fi o ip for some /i G X(S'^) we get 
T^E{fj, oip) = Gdo^ii'Logdo^i c, {d o ipiYu) . 
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After introducing an inner product on X{S^), we could 
compute the gradient of E{tp) by solving 

{V^E, u)x(si) = Grfo¥.i(Logrfo>pi c, {d o ipiYu) . 

There is however a better expression for the gradi- 
ent, which fits better with the action of DifF(S'^) on 

Consider the related function 

E{(p) = i dist(c o (pi, d)"^ . 

Its gradient can be computed using the rnvariance of the 
metric under ip, i.e. E{ip) = E{ip\), which implies 

T^E{Sip) = T^iE {-{ipiYdifi o ipi) 



-T^iE ( ( —Sf I o ipi 



Gdoip (- Logrfoip c, {d! o Lp) 5ip) . 



Using the invariance of Log and the metric Gc under 
reparametrizations, 

Logrfoy c o ^ = (Logrf c)op 
Gdocpih o ip^koip) ^ Gd{h, k) , 

we get 

T^E{6ip) = Gdoif, (- (Log^(c o ip\)) o (d! o -p) 5p) 
= Gd {- Log^(c o pi), d'{5p o pi)) , 

and hence the gradient can be obtained by solving 

{V^E, fi)x(s^) = -Gd(Log^{c o (pi), d'^i) . 

Since the fimctions E and E differ only by exchanging 
c and d, we have another way to express the gradient of 
E, 

{V^E, ti)x{s^) = -Gc(Log^(d o pi),c'fi) . (10) 

As the inner product on the space of vector fields 
we can use the one induced by the identification of 
X(5^) with the vertical space at the curve c using the 
infinitesimal action Cm(c) = c'/i. The inner product is 
given by 

i^>x(si) = Gc(cV, cV) . (11) 
Theorem 6.3. The (right-trivialized) gradient of the energy 

E{p) = ^ dist(c, d o (pi)'^ 

with respect to the inner product | (Tl) is given by the vec- 
tor field on corresponding to the vertical projection of 
— Logg(d o p\), i.e. 

c'W^E = ~ Verc (Log^((i o (pi)) . 

It can be computed by solving the equation 

Gcic'V^E, c'fji) = -Gc(Log,(d o pi),c'^i) 

for V^E e X{S^) with ^ e X{S^) being an arbitrary vector 
field. 



Proof: Note that the equation 

Gc{h, c'fi) = -Gc(Log^ d o (pi, c'fj.) 
is satisfied for all fi G X(S'^) if and only if 

h = - Veic (Log^{d o p^i)) . 

Each element in the vertical space corresponds via the 
infinitesimal action to one element of X{S^). The theorem 
now follows from combining equations | (T0) and | (TT} . □ 
The next algorithm computes the element ii E X{S^) 
corresponding to the orthogonal projection of a vector 
h £ Tcluim{S^ ,R^) to the vertical subspace Ver(c). 
function Ver(c, h) 

V test function on X(S'^) 
[I solution of Gc{d v) = Gc{h,c!v) 
return [i 
end function 

We use a finite element method with Lagrange elements 
of first order to numerically compute the vertical pro- 
jection. A more explicit formula for the inner product 

Gc(ft., fc) is given by 



G,(h,k) 



{h',k') 3{h',c'){k',c') 



d0. 



The algorithm to find geodesies between 
imparametrized shapes C, D takes as input two 
parametrizations c, d of these shapes such that 
7r(c) = G and 7r(c?) = D and finds the diffeomorphism 
lp € Diff (5^), such that the geodesic distance dist(c, doip) 
is minimal. 

function SolveBVP(c, d, N) 

ip Idgi Notation tp := tpi in ||9| 

while dist(c, do tp) is not minimal do 
h ^LOG{c,doip,N) 

^ <-Ver(c, h) > Notation fi -V^E 

rj <-FLOW(-//,Idsi,a) 

end while 
return tp 
end function 

To understand the algorithm note that tp corresponds to 
(pi in Theorem 6.3 There is no need to compute (p itself, 
as only pii is neccessary to compute the reparametriza- 
tion of the curve d. In each iteration of the algorithm we 
first compute the gradient ji = —V^E with the help of 
Theorem 6.3 A continuous gradient descent would take 
the form 

dt(p ^ -VipE o (p . 
A first order time-discretization would be 



FP(a,Idsi)o(^, 



with F1'^(q:, Id5i) denoting the flow of the vector field 
ji — —V^E up to time a starting from Idgi at time 0. 
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Since we are only interested in ip — ipi, we can rewrite 
this as 

ipi+i = ipi+ii = v3iioFl''(a,Idsi)i 
= Vz°Fr^(a,Idsi). 

The gradient descent step is repeated until the relative 
decrease of dist(c, doip) falls below a prescribed thresh- 
old. 

Similarly to the function Logg^^{qi) in Section 



6.1 



the 

parameter a is the step-size for the gradient-descent and 
is chosen adaptively to ensure the distance decreasing in 
each step. 

6.4 Adaptive grid refinement 

The behavior of shortest paths between unparametrized 
curves can be understood as a combination of stretching 
and bending. Bending of a parametrized curve is nu- 
merically well behaved, since uniformly sampled curve 
will stay approximately uniformly sampled. Stretching 
however tends to expand a very short section of the 
curve into a much larger one and will lead to a curve 
that is under-sampled in the expanded area, unless 
the original parametrization was chosen to counter this 
effect. See Figure. |2] for an example of this. 



In the gradient descent algorithm from Section 6.2 



we 



start with a curve that is sampled uniformly at the points 



2tt 



■ 5 — 1 



-2n. 



At each iteration, given a grid xo,...,a;„-i for the 
curve c, the curve d would be sampled at the points 
^{xo), ■ ■ ■ , ^p{xn^i). Stretching of the curve c would cor- 
respond to a large derivative ^"'(2;) or a large distance 
1^(2^4+1) — ij{xi)\ between two consecutive points. We 
add points to the grid, whenever the distance exceeds 
that of a uniform grid, i.e. 



|V'(a;i+i) - i^{xi)\ > 



2n 

no 



Here uq is the size of the original grid. By this we ensure 
that the target curve dot/j will not be under-sampled. 

Similarly we remove a point Xi, whenever the neigh- 
bouring points would be sufficient to provide enough 
resolution, i.e. the conditions 



27r 

\x^+i < — and \ip{xi+i) 

no 

are satisfied. 



4>{xt-i)\ < 



27r 
no 



7 Results 

In this section we present experiments demonstrating 
our numerical computations. The algorithms described 
in Section |6] were implemented in Python using the 
NumPy and SciPy libraries. To compute the vertical 
projection of a tangent vector as described in Section |6.2| 
we used the finite element library FEniCS |13|. The 
shapes used for the experiments come from the database 



of closed binary shapes collected by the LEMS Vision 
Group at Brown university (http: / / www.lems .brown." 
edu/~dmc|. Each curve was initially parametrized by a 
set of 300 points positioned at equal arc length along 
the curve. Each geodesic between curves with fixed 
parametrizations was computed using 25 time steps. To 
find the optimal parametrization the gradient descent 
algorithm from Section 6.2 with a maximum of 100 
iterations was used. 

In Table [T| we show the geodesic distance between 
several shapes, the number of iterations neccessary to 
find the optimal reparametrization and the number of 
points on the discretized curve. Figure [l] shows examples 
of some of the geodesies. The geodesies are sampled 
at timesteps 0, 5, 10, 15, 20 and 25. The first and last 
images in each row thus show the template and target 
curves respectively. While the template curve is always 
parametrized proportional to arc length, the target curve 
will in general have an irregular parametrization, while 
still representing the same shape. 

From a mathematical point of view we would expect 
the distance between two shapes to be symmetric, i.e. 
d(Co, Ci) — d{Ci,Co)- However, if we look at the table, 
we see that the distance from the cat to the dog is 
about 20% smaller than the distance from the dog to 
the cat. The difference is even more striking in the 
distances between the cow and the dog. The reason for 
this difference can be foimd by looking at the minimizing 
geodesies. The optimal geodesic from the dog to the 
cow is shrinking a leg of the dog and at the same time 
growing another leg of the cow from the body. As we 
can see in Figure |3] such growth tends to originate from a 
single point. Thus one point of the dog will be expanding 
to create the whole leg of the cow, while the whole leg 
of the dog will collapse to one point on the cow. This 
is why it is neccessary to dynamically add points to 
the discretization of the template and target curves as 
described in Section l6!4l 

In Figure |2] we show a more pronounced example of 
the pathological behaviour of the metric. We compute 
the geodesic between the two shapes in the first row of 
Figure |2] On the left hand side we can see the minimal 
geodesic computed using the grid refinement method 
described in Section 6.4 The lines connecting the two 
shapes show the paths of the marked points along the 
geodesic. We see that the whole fold of the target shape 
is being created out of a single point of the template 
shape. For comparison we also tried to solve the shape 
matching problem without grid refinement. The result 
can be seen in the right image of the bottom row. The 
template shape is lacking the resolution to resolve the 
fold and therefore the matching is rather unsatisfactory. 

8 Conclusions 

Riemannian metrics on shape space of plane curves 
are of great interest in a wide variety of applications 
in image analysis and computational anatomy. Due to 
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Fig. 1 . Examples of geodesies between several shapes. The geodesies shown are from the eat to the dog in the first 
row, from the dog to the cow in the second, from the cow to the donkey in the third and from the donkey back to the 
cat in the last row thus forming a geodesic quadrangle. The marked points on each geodesic are the paths of initially 
equidistantly spaced points of the template curve showing the stretching and compression along minimal geodesies. 

TABLE 1 

Geodesic distances between several shapes as computed by the minimization algorithm. The number of iterations 
needed to achieve the minimum and the number of points on the refined grid are also shown. From a numerical point 
of view the distinction between template and target is significant. On average the distances differ by about 14%. 
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the infinite dimensional nature of shape space, metrics 
that allow for efficient computations of geodesies and 
distance between shapes are particularly useful. 

In this paper we generalize the i?-transform, first 
introduced in the work of [23 1, to the class of all elastic 
metrics, whose coefficients satisfy 4b^ > a^, where a and 
b are parameters controlling the degree of bending and 
stretching of the curve respectively (see Section [LTJ. This 
transformation allows us to obtain efficient algorithms 
for computing the geodesic distance between shapes as 
well as to gain a better understanding of the geometry of 
the space. For the case of open curves we obtain explicit 
formulas for geodesies and we can show that the space 



of parametrized open curves is a flat space in the sense 
of Riemannian geometry. As a consequence, it follows 
that the shape space of unparameterized open curves has 
positive sectional curvature. For closed curves the situ- 
ation is more difficult, since the space of parametrized 
curves is not flat anymore, but the representation is still 
very useful from a numerical point of view. 

We have presented experiments showing the minimal 
geodesies between closed, unparametrized curves and 
computed the geodesic distances between them. For 
certain pairs of shapes the i?-transform allows for very 
efficient computations, while for other pairs we see that 
the tendency of the metric to compress and stretch points 
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Fig. 3. The top rows shows the minimal geodesic from the shark and the plane. The bottom row shows the minimal 
geodesic computed with template and target curves reversed. We see that the intermediate shapes almost coincide. 




Fig. 4. Similarly as in Figure|3]the top row shows the minimal geodesic from the cat to the dog. The bottom row shows 
the minimal geodesic computed with template and target curves reversed. Since in this case there is stretching and 
compression happening along the geodesies the two paths are not symmetric any more. This also corresponds to the 
difference in the geodesic distances as can be seen in Table[l] 



along the curve leads to numerical difficulties. In particu- 
lar the difference in the geodesic distance, depending on 
which curve was chosen as template and which as target 
shows that more work is required to understand this 
behaviour and to develop numerical methods capable 
of dealing with this situation. We plan to address these 
questions in future work. 

Another way to address this problem would be to 
use metrics of higher order, i.e. those that depend on 
higher derivatives of either the curve c or the variation h. 
Finding efficient representations of such metrics similar 
to the i?-transform used for metrics of order one remains 
an open question. 
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